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P^; ABSTRACT 

'^N \ Rigorous electromagnetic field simulations are an essential part for scatterometry and mask pattern design. 
\^ • Today mainly periodic structures are considered in simulations. Non-periodic structures are typically modeled 
by large, artificially periodified computational domains. For systems with a large radius of influence this leads to 
very large computational domains to keep the error sufficiently small. In this paper we review recent advances 
C/3 ] in the rigorous simulation of isolated structures embedded into a surrounding media. We especially address the 
. Sh ' situation of a layered surrounding media (mask or wafer) with additional infinite inhomogeneities such as resist 
"^ I lines. Further we detail how to extract the far field information needed for the aerial image computation in the 
r-s ' non-periodic setting. 

c/3 ' Keywords: finite elements, transparent boundary conditions, Wiener-Hopf integral. Maxwell's equation, scat- 
. rH ' terometry, microlithography 

>^' 
^ ■ 1. INTRODUCTION 

Cl(i . . . 

The simulation setting we focus on is as depicted in Figure [TJ An incident field Einc enters the computational 
domain from above and excites a scattered field Egc. For simplicity we assume that the incident field is a 
►^ , plane wave. A more complex illumination can be modeled by a decomposition of the source field into plane 
r^^, ' waves. Afterwards the scattered fields of the various plane waves are superimposed accordingly to their mutual 
00 , coherence relation. In Figure [T] the scatterer is embedded into a layered media and completely contained in the 
•/^ ' computational domain. Using a non-periodic domain allows for the reduction of the computational domain to 
^^ , the minimum dimensions of the considered structure. Furthermore, no cutoff error is introduced in contrast 
^^ ' to modeling the structure as a large, artificially periodified computational domain. However, the numerics is 
C3 I more delicate for non-periodic computational domains, since special transparent boundary conditions must be 
implemented. Besides the more involved numerics for the near field computation the extraction method of the 
far field data must also account for the non-periodic setting. In a non-periodic direction the diffraction modes 
^ ' are not discrete anymore. Instead we must deal with a band-limited continuous Fourier transform. 

J^ ' In earlier papers^' ^ we already studied the situation as in Figure [T] In this paper we address the more com- 

H I plicated situation as in Figure [21 Now, the scatterer is not completely enclosed by the computational domain 
but is itself infinitely prolonged as indicated by the dashed lines. Such configurations are called "waveguide 
inhomogeneities".^''* The theory presented in our papers'^'"* assumes that the incoming field is zero on the 
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Figure 1. Sketch of the simulation setting. An incident field Einc enters the computational domain from above and excites 
a scattered field Eac- The computational domain encloses the scatterer. The surrounding media is typically layered. More 
general situations are depicted in Figure [2] 




Figure 2. More complicated geometries. We aim to compute the field scattered off an infinite resist line's tip (left hand 
side) or off an infinite plate's corner. Periodicity in one coordinate direction may also be incorporated but is not especially 
treated in this paper. 

infinite waveguides outside the computational domain. Therefore, deahng with incident plane waves necessi- 
tates a non-trivial extension of the scattering problem formulation which is presented in the next section. The 
presented theory generaUzes classical results by Wiener and Hopf for ideal geometries like a thin half-screen, 
see. ^'^Accordingly, a new far field extraction formula must be derived for this setting. Numerical aspects on this 
novel approach are discussed in Section [31 



2. ELECTROMAGNETIC SCATTERING PROBLEMS 

2.1 Maxwell's equations 

Starting from Maxwell's equations in a medium without sources and free currents and assuming time-harmonic 
dependence with angular frequency uj > the electric and magnetic fields 



must satisfy 



E(a;, y, z, t) = E(a;, y, z)e "^■*, H(x, y, z, t) = H(x, y, z)e'~ 



V X E = iw/iH, V ■ eE = 0, 

V X H = -icjeE, V ■ ^iH = 0. 



Here e denotes the permittivity tensor and /j, denotes the permeabiHty tensor of the materials. In the following we 
drop the wiggles, so that E ^ E, H -^ H. From the equations above we then may derive (by direct substitution) 
the second order equations for the electric field 

V X /i^^V X E - Lu^eE = 0, 
V • eE = 0. 

Similar equations hold true for the magnetic field - one only need to replace E by H and interchange e and 
H- Observe that any solution to the first equation also meets the divergence condition (second equation). We 
therefore drop the second equation in the following. 
For the sake of a simpler and more compact notation we rewrite these equations in differential form, 

din^^die- uj'^ee == 0. (1) 

A reader not familiar with this calculus may replace the exterior derivatives do, di, d2 with classical differential 
operators, do —> V, di — > Vx and ^2 — > V-. Here, the electric field appears as a differential 1-form, e = 
Cxdx + Cydy + Czdz, whereas the material tensors act - from a more mathematical point of view - as operators 
from the differential I-forms to the differential 2- forms, 

e, /i : Alt^ -^ Ah^. 

In the following we drop the sub- indices for the exterior derivatives do , di , ^2 ■ 

2.2 Simple structured exterior domains 

We consider a situation of an isolated scatterer as in Figure [TJ In a scattering problem one compares a reference 
configuration with a perturbed configuration containing the scatterer. In Figure [3] the reference configuration is 
a wafer/mask blank with its layered media. For a prescribed illuminating plane wave the solution to Maxwell's 
equations in the layered media is denoted by Erof. Hence the field E^ef already comprises the reflections at the 
layered media. The point here is that Ercf is a solution to Maxwell's equations for the reference configuration and 
can be directly computed. Now assume that the computational domain contains a scatterer as in Figure [3] (right). 
Actually the computational domain may contain any change of the material properties. For the prescribed 
illuminating plane wave we denote the total electric field by Etot and introduce the scattered field as 

Esc = Etot — Erof. 

The scattered field Ego is the difference between the reference field and the field observed for the perturbed 
geometry. This way the scattered field Ego contains no incoming field components, so that a radiation boundary 
condition can be imposed. This yields the following formulation of the scattering problem as a coupled interior- 
exterior domain problem, 

d/i~ de — uj ee = in SI 

d/i^ deso — w eesc = in R \ Vt 

e = Gso + Giof on 9f7 

dfx^^e = diJ,~^esc + d/i^^Giof on dfl. 



Here, fi denotes the computational domain and e the total field Etot in the interior domain. The first equation is 
Maxwell's equation in the interior domain for the total electric field. The second equation is Maxwell's equation 
for the scattered electric field in the exterior domain. The last two equations are the coupling conditions on 
the boundary. The third equation enforces tangential continuity of the total electric field, whereas the fourth 
equation enforces tangential continuity of the magnetic field. As mentioned above a radiation boundary condition 
must be imposed for the scattered field Sgc- We postpone the formulation of the radiation boundary condition 
to Section [21 
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Figure 3. Scattering problem. For an illuminating plane wave, Maxwell's equations are solved for the reference configura- 
tion (left figure) to obtain E,of ■ The perturbed geometry with the scatterer is shown on the right hand side. The scattered 
field Esc is the diff'erence between the total field Etot and the reference field, Esc = Etot — E,.cf . The scattered field satisfies 
a radiation boundary condition. For the formulation of the scattering problem the reference field Ej-cf is needed on the 
boundary of the computational domain. 
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Figure 4. 2D scattering problem with complicated exterior domain. The exterior domain has different layered media on 
the left and right side. One might use the reference geometry on the left hand side, but solving for the reference field Ercf 
is not trivial in this case. 



2.3 Multiple structured exterior domains 

We now consider the more complicated situation from Figure [2l For the 2D setting we refer to Figure H) The 
exterior domain in Figure [4] is differently layered on the left and right horizontal half space. Hence it is not clear 
how to define the reference geometry. One might arbitrarily use a reference domain as given in Figure 0] (right), 
but this leads to a reference problem with a non-trivial solution. Here we propose a different approach. We 
horizontally split the overall space as in Figure [5] (left) to get two reference solutions Ejcf^i and Ercf, 2 in the 
different layered media for a prescribed illuminating plane wave. On the right hand side of Figure the exterior 
domain is splitted into two parts i~lcxt,i and r2cxt,2- Now we deal with two different scattered fields defined by 



Esc,l 
Esc, 2 



Etot — Ei-cf,i 

Etot — Ei-cf,2, 



(2a) 
(2b) 



where Etot is again the total field. The scattered fields Esc,i and Esc,2 satisfy Maxwell's equations on ricxt.i 
and r2cxt,2 respectively and are outward radiating. From the definitions ^ we conclude that the scattered fields 
jump across the interface F = ilcxt,i H 51cxt,2, 



dfi^'^'Es^ 



Esc, 2 - 

I -d^i' 



lE 



sc,l — Ercf,! — Ei.cf,2 

sc,l — dfl Ei-of,! — dfj, Ercf, 2- 



Note that this jump is exactly the difference of the reference fields. Since for both reference fields the same 
incoming plane wave is chosen the difference is a purely outgoing wave. This fact will play an important role, 




J^n 



^E,e,l^ 



^ Esc,2 4< 
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Figure 5. Choice of different reference fields for the splitted exterior domains. Ercf,i and Erof,2 are solutions to Maxwell's 
equation in the different layered media for a prescribed illuminating plane wave. On the right hand side the exterior domain 
is splitted into two parts ficxt.i and ilcxt,2. On each of the domain we define scattered fields Esc.i, Esc, 2 as the difference 
of the total field Etot and Erof.i, Erof,2 respectively. Additional coupling conditions on the boundary F — Jlcxt.i n ficxt,2 
are needed. 



when we define transparent boundary conditions. The scattering problem with two different reference fields now 
reads as, 



dfi de — uj ee 

dfi~^desc,k - w^eescfe 

e 

d/j^^e 

*^sc.2 6sc,l 



inf7 

in l^cxt.k 

esc.fe + ercf ,fc on dil n dflcxt,k 

d^^^Bsck + dfi^^e,:ct,k on dfl n 917cxt,k 

Grcf,! — ei.cf,2 on 951oxta n 9Sloxt,2 

d^^^esc,2 - dfl^^ Bscl = d^^^Grcf,! - rf/i"^ercf,2 On 9Sloxta H 9rioxt,2 



(3a) 
(3b) 
(3c) 
(3d) 
(3e) 
(3f) 



for k ~ 1,2. The generalization to muhiple reference solutions is straightforward. One only needs to account 
for all coupling conditions at the interfaces of two adjacent exterior domains. As a strategy one decompose the 
exterior domain such that the reference fields to the prescribed incoming field are numerically available. In the 
case of a layered media, the reference solution can be obtained quasi-analytically. For waveguide inhomogeneities 
as in Figure [5] (left) one may also split the domain at y = 0. The domain with y < is a simple layered media. 
The domain with y > contains the waveguide (resist line) and the corresponding reference solution can be 
obtained by solving a 2D scattering problem on the cross-section y — const. However, there is a restriction for 
this construction. The difference of two reference fields in the coupling condition ^ must satisfies a radiation 
boundary condition defined in the next section. 

2.4 Weak formulation 

The scattering problem ^ is given in strong form. To use the finite element method one must cast the scattering 
problem into weak form. Since the problem is still posed on the entire space R'^ we use test functions with local 
support. In the following we apply the partial integration rule for differential 1-forms u, v, 



u A rfv = / du A V 
n Jn 



u A V. 



So given a test function v we get. 



dv /\ 1.1 de — oj V A ee 



on 



an 



V A /i ^de 



/ dv/i ^descfc - tj^eescfc = / v A /i ^desck, k ^ 1,2. 



Summing these three equations and using equations (pd]) . jSj yields 

/ dv A ^^"^de — tj^v A ee = / v A ^^"^deref.i + / v A /x^"^deref,2 + 

/ V A^"^d(e,.of,2 -erof,i) , 

v'aOcxt,inosicxt,2 

where e denotes the compound field e = e + egci +esc,2- This equation accounts for the coupling of the Neumann 
data on the various boundaries. But e jumps across the boundaries, so that it is not feasible in a finite element 
discretization. In a further step we therefore must incorporate the Dirichlet coupling conditions (|3d[) and ((3^ . 
To do that one adds supplemental fields with local support near the interface boundaries to e making the so 
constructed field e' continuous across the boundaries. The choice of the additional fields is not unique, but 
must account for the construction of transparent boundary conditions in the next section. The corner point 
dQ n (?r2oxt,i n rioxt,2 needs a special treatment for a proper matching of the Dirichlet values at this point. Let us 
denote by i?e a small vicinity of this corner point. We define a smooth function x with the following properties, 

X = 1 on dil 

X ^ on9fiext,i n5rJcxt,2\S<:. 

Now we introduce fields ej.^j -^ and e'^^^ ^ such that 

Kel,2 = Xercf,2 + (1 -x)(erof,2 -erof,l). 

This allows for the definition 

e' = e + e;.^f i+e;„f 2, 

which is continuous across the coupling boundaries. The variational problem for e' now reads as 



/ dv A fi ^de' — uj'^v A ee' = / v A /i ^derci.i + v A /i ^derci.2 + 

/ V A/^"^d(eref,2 -ercf,i) + 

Jc5ncxt,inc?ncxt,2 
/ dv A ^"-^(iej.(,f 1 - cj^v A eej.gf 1 + 

/ dvA/^-ide;.,f2-'^^vAee;,^f2. (4) 

The construction of e' looks very tricky in this analytic setting, but is very simple in the finite element context. 
An explicite construction of the function x can be avoided. Instead one imposes the Dirichlet jump conditions 
directly on the finite element coefhcient vector. 

3. TRANSPARENT BOUNDARY CONDITIONS/PERFECTLY MATCHED LAYERS 

So far, in all considerations the various scattering problems were posed on the entire domain R'^ and are therefore 
numerically not feasible. This is overcome by using transparent boundary conditions. We use the perfectly 
matched layer method introduced by Berenger.'''^ This method exploits the analytic continuation properties of 
the scattered field in the exterior domain. In a nutshell using an appropriate complex continuation, the scattered 
field is transformed to an exponentially decaying field without affecting the matching condition with the field in 
the interior domain. This allows for a truncation of the computational domain. In an earlier paper we proposed 
an extremely efficient adaptive PML method which also copes with various kinds of inhomogeneous exterior 
domains. ^^^ The method can also be applied to Q. One derives a weak form of Maxwell's equation with PML 
which has exactly the same form as equation Q, but where the fields erof,i, ercf,2, Ggci and esc,i are replaced by 



led = 310nm 



zk 




hmosi = 65nm 



Figure 6. Test problem with different right and left hand side exterior domains. A MoSi structure (nmosi = 2.52 + 0.596i) 
is mounted on a glass substrate (wg = 1.5306). The initial finite element mesh is shown. Quadrilaterals are modeled 
as infinitely prolonged by using the adaptive PML method. The structure is illuminated by a s-polarized plane wave 
{Ex, Ez — 0) with vaccum wavelength A = 193nm and wave direction k = ( — 1/\/2, 0, — 1/%/2)- To compare the simulated 
field with the periodic setting in Figure [7] the field values are extracted on the dashed lines. 
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Figure 7. Same setting as in Figure [S] but periodic boundary conditions are used in horizontal direction. The periodicity 
length Itt was gradually increased to suppress the cut-off error caused by the artificial periodification. Field values are 
extracted on the dashed lines placed as in Figure [6] 

their complex continuation. The last four terms on the right hand side of equation Q are non-standard. These 
terms model the coupling of the illuminating field with non-trivial structures in the exterior domain. Within the 
PML the magnitude of these further right hand side terms is proportional to ||ei.ef.2 — ei.ef,i||. As far as eicf,i and 
ercf,2 share the same incoming field, we conclude that the additional terms are also evanescent within the PML. 
This must be guaranteed when splitting the exterior domain into different reference geometries. 



4. IMAGING 

In this section we briefly sketch how to extract the far field information. We first refer to Figure [3] with of 
a simple layered media in the exterior domain. A rigorous field evaluation formula of the scattered field egc 
in the upper half space can be derived by means of the Rayleigh-Sommerfeld integral. ^'^ In microlithography 
application one is mainly interested in the far field data. Moving an evaluation point towards infinity in a 
direction r(nx, ny, nz), t -^ cxd, one shows that the field value converges up to a scaling to the Fourier transform 
of the scattered field in the xy— plane. 

Alternatively to an evaluation of the Rayleigh-Sommerfeld integral, one may use Green's tensor Gj , 



dfJ. ^dgSp 



- w^egs^ 



S{. 



in R3 \ 17 



in the exterior domain to evaluated the field in the exterior domain. One gets 



z(Xp) 



cM ^^Xjj Sajp/^ Q^Ggc, 



an 



so that the scattered field is only needed on the boundary of the computational domain. Since typically one is 
only interested in the far field, it is sufficient to use the asymptotics of the Green function for Xp = T{nx, ny, nz), 
T — * CXD, which can be easily computed for the layered media. The situation for the more complicated exterior 
domain in Figure [His more involved. The computed scattered field Ggc depends on the choice of reference field 
Ercf • The reference field already comprises refiected fields from the infinite resist line. From a metrology point of 
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Figure 8. Convergence of the periodic solution with increasing periodicity length l-n to the solution obtained with the new 
method. The field values were extracted on horizontal lines above and below the structure as in Figures [H] [T] Even for 
moderate accuracy demands one observes a tremendous speed-up by the new method compared to the periodic setting. 

view this corresponds to a calibration. Once a reference field in the entire exterior domain R^ \ J7 as in Figure 3] 
(left) is fixed, the scattered field egc = etot — ^rci has a purely continuous spectrum and can be computed by 
means of either a Rayleigh-Sommerfeld integral or by using of Green's tensor asymptotics. The computation of 
the Green function's asymptotics requires a further numerical simulation of a scattering problem for each far 
field evaluation point. We detail this in a later paper. 

5. NUMERICAL EXAMPLE 

We applied the new method to the simple 2D test problem given in Figure [S] An analytic solution to this 
problem is not known. To verify the accuracy of the method we simulated the structure on an artificially 
periodified domain as in Figure El When increasing the periodicity length Z^ one expects convergence to the 
solution obtained for the isolated case. To measure the error we extracted the electric field values on horizontal 
lines below and above the structure (Figures [SI [7]) • As can be seen from Figure [5] the field values for the periodic 
setting converges to the values obtained with the new method confirming our theory of scattering within multiple 
exterior domains. Even for moderate accuracy demands one observes a tremendous speed-up by the new method 
compared to the periodic setting. 
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